setwd("C:/Users/vinic/Dropbox/ABCP2020")
library(rio)
data <- import(file="Dados_ABCP2020.xlsx")


#Construindo as vari�veis

#engajamento de alta intensidade (dependente)
engajamento <- data$P21
engajamento <- as.factor(engajamento)
table(engajamento)
levels (engajamento) <- c("0","0","0","0","0","0","1")

#escolaridade (1 para universit�rio completo/incompleto)
escolaridade <- data$P5
escolaridade <- as.factor(escolaridade)
table(escolaridade)
length((escolaridade))
levels(escolaridade) <- c("0","0","0","0","0","0","0","1","1")

#sexo (1 para homem)
sexo <- data$P13
sexo <- as.factor(sexo)
table(sexo)
levels(sexo) <- c("0","1","0","0")

#ideologia do filiado (centro-direita como 1 - ponto 6-10 da escala)
ideologia <- data$P50
ideologia <- as.factor(ideologia)
table(ideologia)
length(ideologia)
levels(ideologia) <- c("0","0","0","0","1","1",
                       "1","1","0","1","0","0","0")


#motivo de filia��o - oposi��o (para se opor a um partido como 1 - muito importante motivo para filia��o)
motivo_oposicao <- data$motivo_opor
motivo_oposicao <- as.factor(motivo_oposicao)
table(motivo_oposicao)
length(motivo_oposicao)
levels(motivo_oposicao) <- c("0","1","0","0")


#cargo no partido (militante apenas, codificado como 1)
cargo_part <- data$P26
cargo_part <- as.factor(cargo_part)
table(cargo_part)

length(cargo_part)
levels(cargo_part) <- c ("0","0","0","0","0","0","0",
                         "0","0","0","0","1","0","0",
                         "0","0","0","0","0","0","0",
                         "0","0","0","0")



#### modelo #### 
modelo1 <- glm(engajamento ~ sexo + escolaridade + ideologia +
                 cargo_part + motivo_oposicao, family="binomial")
summary(modelo1)


#### interpretando os coeficientes ####

#motivo_oposicao
exp(0.24494) #extraindo a raz�o de chance
#1.277545

#sexo
exp(0.39947)
#1.491034

exp(0.02843)
#1.028838 (obs.: escolaridade n�o � significativa)

#ideologia
exp(-0.29194)
#0.7468133

#apenas militante (sem cargo no partido)
exp(-0.47886)
rc <-0.6194892
p_semcargo <- rc/(1+rc)
p_semcargo

p_cargo <- 1-p_semcargo
p_cargo

#probabilidade de ter um cargo no partido e ser altamente engajado � de 0.6174786,
#enquanto a de n�o ter se reduz a 0.3825214 - raz�o de chance em favor de ter cargo (1.614233)
r <- 0.6174786/0.3825214
r
#1.614233



#### EXPERIMENTO SUJEITO �NICO ####

#An�lise do experimento considerando apenas quem participou de todos os cen�rios #1917 respondentes) 

dados <- import(file="dados-unesp-limpo.xlsx")

#cen�rio controle
eng_contr <- dados$hs_controle_cat
eng_contr <- as.factor(eng_contr)
table(eng_contr)

#cen�rio vit�ria do pr�prio partido
eng_ben <- dados$hs_benef_cat
eng_ben <- as.factor(eng_ben)

#cen�rio vit�ria do partido que odeia
eng_odio <- dados$hs_odio_cat
eng_odio <- as.factor(eng_odio)

#cen�rio democracia interna
eng_demint <- dados$hs_demint_cat
eng_demint <- as.factor(eng_demint)

#descri��o do engajamento de alta intensidade nos cen�rios controle, vit�ria do pr�prio partido,
# �dio e democracia interna
length(eng_contr) #n=1917
table(eng_contr=="7") # alta intensidade controle (366/1917) = 0.1909233
table(eng_ben=="7") # alta intensidade vit�ria do pr�prio partido (800/1917) = 0.4173187
table(eng_odio=="7") # alta intensidade �dio (702/1917) = 0.3661972
table(eng_demint=="7") # alta intensidade democracia interna (615/1917) = 0.3208138

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr=="7", eng_ben=="7", alternative = "less", paired = T) #p-value < 2.2e-16

#efeito do engajamento por �dio
t.test(eng_contr=="7", eng_odio=="7", alternative = "less", paired = T) #p-value < 2.2e-16

#efeito do cen�rio democracia interna
t.test(eng_contr=="7", eng_demint=="7", alternative = "less", paired = T) #p-value < 2.2e-16

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
#(o primeiro tem intensidade maior?)
t.test(eng_ben=="7", eng_odio=="7", alternative = "greater", paired = T) #p-value = 6.808e-08***


####comparando engajados �dio e antipetismo (base apenas participantes do experimento de sujeito �nico) ####
altamente_odio <- dados$hs_odio
altamente_odio <- as.factor(altamente_odio)
table(altamente_odio)
length(altamente_odio)
levels(altamente_odio) <- c("0","0","0","0",
                            "0","0","1")


antiPT <- dados$rej_PT
antiPT <- as.factor(antiPT)
table(antiPT)
length(antiPT)
levels(antiPT) <- c("0","1")


table(antiPT, altamente_odio)
#fen�meno mais amplo do que antipetismo (dos altamente engajados, apenas 35,04% rejeita PT: 246/702)


#### subsets para examinar resultados do exp de sujeito �nico por partido ####

#### PT ####
dadosPT <- dados[which(dados$partido=="PT"),]
eng_contr_PT <- dadosPT$hs_controle_cat
eng_contr_PT <- as.factor(eng_contr_PT)

eng_ben_PT <- dadosPT$hs_benef_cat
eng_ben_PT <- as.factor(eng_ben_PT)

eng_odio_PT <- dadosPT$hs_odio_cat
eng_odio_PT <- as.factor(eng_odio_PT)

eng_demint_PT <- dadosPT$hs_demint_cat
eng_demint_PT <- as.factor(eng_demint_PT)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PT) #n=188
table(eng_contr_PT=="7") # alta intensidade controle (42/188) = 0.2234043
table(eng_ben_PT=="7") # alta intensidade benef�cio pr�prio partido (89/188) = 0.4734043
table(eng_odio_PT=="7") # alta intensidade �dio (76/188) = 0.4042553
table(eng_demint_PT=="7") # alta intensidade democracia interna (60/188) = 0.3191489

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PT=="7", eng_ben_PT=="7", alternative = "less", paired = T) #p-value = 4.038e-12

#efeito do engajamento por �dio
t.test(eng_contr_PT=="7", eng_odio_PT=="7", alternative = "less", paired = T) #p-value = 2.136e-07

#efeito do demint
t.test(eng_contr_PT=="7", eng_demint_PT=="7", alternative = "less", paired = T) #p-value = 0.00163

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PT=="7", eng_odio_PT=="7", alternative = "greater", paired = T) #p-value = 0.004483 ***


#### PSDB ####
dadosPSDB <- dados[which(dados$partido=="PSDB"),]
eng_contr_PSDB <- dadosPSDB$hs_controle_cat
eng_contr_PSDB <- as.factor(eng_contr_PSDB)

eng_ben_PSDB <- dadosPSDB$hs_benef_cat
eng_ben_PSDB <- as.factor(eng_ben_PSDB)

eng_odio_PSDB <- dadosPSDB$hs_odio_cat
eng_odio_PSDB <- as.factor(eng_odio_PSDB)

eng_demint_PSDB <- dadosPSDB$hs_demint_cat
eng_demint_PSDB <- as.factor(eng_demint_PSDB)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PSDB) #n=125
table(eng_contr_PSDB=="7") # alta intensidade controle (17/125) = 0.136
table(eng_ben_PSDB=="7") # alta intensidade ben proprio part (39/125) = 0.312
table(eng_odio_PSDB=="7") # alta intensidade odio (39/125) = 0.312
table(eng_demint_PSDB=="7") # alta intensidade demint (33/125) = 0.264

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PSDB=="7", eng_ben_PSDB=="7", alternative = "less", paired = T) #p-value = 3.289e-05

#efeito do engajamento por �dio
t.test(eng_contr_PSDB=="7", eng_odio_PSDB=="7", alternative = "less", paired = T) #p-value = 5.571e-05

#efeito do demint
t.test(eng_contr_PSDB=="7", eng_demint_PSDB=="7", alternative = "less", paired = T) #p-value = 0.001096

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio
t.test(eng_ben_PSDB=="7", eng_odio_PSDB=="7", alternative = "greater", paired = T) #p-value = 0.5 !(m�dias iguais)




#### MDB ####
dadosMDB <- dados[which(dados$partido=="MDB"),]
eng_contr_MDB <- dadosMDB$hs_controle_cat
eng_contr_MDB <- as.factor(eng_contr_MDB)

eng_ben_MDB <- dadosMDB$hs_benef_cat
eng_ben_MDB <- as.factor(eng_ben_MDB)

eng_odio_MDB <- dadosMDB$hs_odio_cat
eng_odio_MDB <- as.factor(eng_odio_MDB)

eng_demint_MDB <- dadosMDB$hs_demint_cat
eng_demint_MDB <- as.factor(eng_demint_MDB)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_MDB) #n=161
table(eng_contr_MDB=="7") # alta intensidade controle (28/161) = 0.173913
table(eng_ben_MDB=="7") # alta intensidade ben proprio part (53/161) = 0.3291925
table(eng_odio_MDB=="7") # alta intensidade odio (49/161) = 0.3043478
table(eng_demint_MDB=="7") # alta intensidade demint (43/161) = 0.2670807

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_MDB=="7", eng_ben_MDB=="7", alternative = "less", paired = T) #p-value = 7.43e-05

#efeito do engajamento por �dio
t.test(eng_contr_MDB=="7", eng_odio_MDB=="7", alternative = "less", paired = T) #p-value = 0.000327

#efeito do demint
t.test(eng_contr_MDB=="7", eng_demint_MDB=="7", alternative = "less", paired = T) #p-value = 0.003339

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio
t.test(eng_ben_MDB=="7", eng_odio_MDB=="7", alternative = "greater", paired = T) #p-value = 0.2335 (n est signf)


#### DEM ####
dadosDEM <- dados[which(dados$partido=="DEM"),]
eng_contr_DEM <- dadosDEM$hs_controle_cat
eng_contr_DEM <- as.factor(eng_contr_DEM)

eng_ben_DEM <- dadosDEM$hs_benef_cat
eng_ben_DEM <- as.factor(eng_ben_DEM)

eng_odio_DEM <- dadosDEM$hs_odio_cat
eng_odio_DEM <- as.factor(eng_odio_DEM)

eng_demint_DEM <- dadosDEM$hs_demint_cat
eng_demint_DEM <- as.factor(eng_demint_DEM)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_DEM) #n=88
table(eng_contr_DEM=="7") # alta intensidade controle (22/88) = 0.25
table(eng_ben_DEM=="7") # alta intensidade ben proprio part (37/88) = 0.4204545
table(eng_odio_DEM=="7") # alta intensidade odio (35/88) = 0.3977273
table(eng_demint_DEM=="7") # alta intensidade demint (33/88) = 0.375

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_DEM=="7", eng_ben_DEM=="7", alternative = "less", paired = T) #p-value = 0.0007044

#efeito do engajamento por �dio
t.test(eng_contr_DEM=="7", eng_odio_DEM=="7", alternative = "less", paired = T) #p-value = 0.001986

#efeito do demint
t.test(eng_contr_DEM=="7", eng_demint_DEM=="7", alternative = "less", paired = T) #p-value = 0.01047

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio
t.test(eng_ben_DEM=="7", eng_odio_DEM=="7", alternative = "greater", paired = T) #p-value = 0.2979 (n est signf)


#### PP ####
dadosPP <- dados[which(dados$partido=="PP"),]
eng_contr_PP <- dadosPP$hs_controle_cat
eng_contr_PP <- as.factor(eng_contr_PP)

eng_ben_PP <- dadosPP$hs_benef_cat
eng_ben_PP <- as.factor(eng_ben_PP)

eng_odio_PP <- dadosPP$hs_odio_cat
eng_odio_PP <- as.factor(eng_odio_PP)

eng_demint_PP <- dadosPP$hs_demint_cat
eng_demint_PP <- as.factor(eng_demint_PP)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PP) #n=106
table(eng_contr_PP=="7") # alta intensidade controle (21/106) = 0.1981132
table(eng_ben_PP=="7") # alta intensidade benef�cio pr�prio partido (29/106) = 0.2735849
table(eng_odio_PP=="7") # alta intensidade �dio (31/106) = 0.2924528
table(eng_demint_PP=="7") # alta intensidade democracia interna (25/106) = 0.2358491

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PP=="7", eng_ben_PP=="7", alternative = "less", paired = T) #p-value = 0.04407

#efeito do engajamento por �dio
t.test(eng_contr_PP=="7", eng_odio_PP=="7", alternative = "less", paired = T) #p-value = 0.02032

#efeito do demint
t.test(eng_contr_PP=="7", eng_demint_PP=="7", alternative = "less", paired = T) #p-value = 0.1868 (n est signf!!)

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PP=="7", eng_odio_PP=="7", alternative = "less", paired = T) #p-value = 0.2977 (n est sign)


#### PL ####
dadosPL <- dados[which(dados$partido=="PL"),]
eng_contr_PL <- dadosPL$hs_controle_cat
eng_contr_PL <- as.factor(eng_contr_PL)

eng_ben_PL <- dadosPL$hs_benef_cat
eng_ben_PL <- as.factor(eng_ben_PL)

eng_odio_PL <- dadosPL$hs_odio_cat
eng_odio_PL <- as.factor(eng_odio_PL)

eng_demint_PL <- dadosPL$hs_demint_cat
eng_demint_PL <- as.factor(eng_demint_PL)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PL) #n=73
table(eng_contr_PL=="7") # alta intensidade controle (14/73) = 0.1917808
table(eng_ben_PL=="7") # alta intensidade benef�cio pr�prio partido (33/73) = 0.4520548
table(eng_odio_PL=="7") # alta intensidade �dio (25/73) = 0.3424658
table(eng_demint_PL=="7") # alta intensidade democracia interna (25/73) = 0.3424658

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PL=="7", eng_ben_PL=="7", alternative = "less", paired = T) #p-value = 5.879e-06

#efeito do engajamento por �dio
t.test(eng_contr_PL=="7", eng_odio_PL=="7", alternative = "less", paired = T) #p-value = 0.007661

#efeito do demint
t.test(eng_contr_PL=="7", eng_demint_PL=="7", alternative = "less", paired = T) #p-value = 0.00338

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PL=="7", eng_odio_PL=="7", alternative = "greater", paired = T) #p-value = 0.01576


#### PSL ####
dadosPSL <- dados[which(dados$partido=="PSL"),]
eng_contr_PSL <- dadosPSL$hs_controle_cat
eng_contr_PSL <- as.factor(eng_contr_PSL)

eng_ben_PSL <- dadosPSL$hs_benef_cat
eng_ben_PSL <- as.factor(eng_ben_PSL)

eng_odio_PSL <- dadosPSL$hs_odio_cat
eng_odio_PSL <- as.factor(eng_odio_PSL)

eng_demint_PSL <- dadosPSL$hs_demint_cat
eng_demint_PSL <- as.factor(eng_demint_PSL)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PSL) #n=60
table(eng_contr_PSL=="7") # alta intensidade controle (8/60) = 0.1333333
table(eng_ben_PSL=="7") # alta intensidade benef�cio pr�prio partido (17/60) = 0.2833333
table(eng_odio_PSL=="7") # alta intensidade �dio (21/60) = 0.35
table(eng_demint_PSL=="7") # alta intensidade democracia interna (20/60) = 0.3333333

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PSL=="7", eng_ben_PSL=="7", alternative = "less", paired = T) #p-value = 0.002822

#efeito do engajamento por �dio
t.test(eng_contr_PSL=="7", eng_odio_PSL=="7", alternative = "less", paired = T) #p-value = 0.001092

#efeito do demint
t.test(eng_contr_PSL=="7", eng_demint_PSL=="7", alternative = "less", paired = T) #p-value = 0.001905

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PSL=="7", eng_odio_PSL=="7", alternative = "less", paired = T) #p-value = 0.1043 (n est signf)


#### PDT ####
dadosPDT <- dados[which(dados$partido=="PDT"),]
eng_contr_PDT <- dadosPDT$hs_controle_cat
eng_contr_PDT <- as.factor(eng_contr_PDT)

eng_ben_PDT <- dadosPDT$hs_benef_cat
eng_ben_PDT <- as.factor(eng_ben_PDT)

eng_odio_PDT <- dadosPDT$hs_odio_cat
eng_odio_PDT <- as.factor(eng_odio_PDT)

eng_demint_PDT <- dadosPDT$hs_demint_cat
eng_demint_PDT <- as.factor(eng_demint_PDT)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PDT) #n=118
table(eng_contr_PDT=="7") # alta intensidade controle (24/118) = 0.2033898
table(eng_ben_PDT=="7") # alta intensidade benef�cio pr�prio partido (59/118) = 0.5
table(eng_odio_PDT=="7") # alta intensidade �dio (51/118) = 0.4322034
table(eng_demint_PDT=="7") # alta intensidade democracia interna (41/118) = 0.3474576

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PDT=="7", eng_ben_PDT=="7", alternative = "less", paired = T) #p-value = 2.731e-09

#efeito do engajamento por �dio
t.test(eng_contr_PDT=="7", eng_odio_PDT=="7", alternative = "less", paired = T) #p-value = 3.688e-06

#efeito do demint
t.test(eng_contr_PDT=="7", eng_demint_PDT=="7", alternative = "less", paired = T) #p-value = 0.002366

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PDT=="7", eng_odio_PDT=="7", alternative = "greater", paired = T) #p-value = 0.05135***




#### PSB ####
dadosPSB <- dados[which(dados$partido=="PSB"),]
eng_contr_PSB <- dadosPSB$hs_controle_cat
eng_contr_PSB <- as.factor(eng_contr_PSB)

eng_ben_PSB <- dadosPSB$hs_benef_cat
eng_ben_PSB <- as.factor(eng_ben_PSB)

eng_odio_PSB <- dadosPSB$hs_odio_cat
eng_odio_PSB <- as.factor(eng_odio_PSB)

eng_demint_PSB <- dadosPSB$hs_demint_cat
eng_demint_PSB <- as.factor(eng_demint_PSB)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PSB) #n=128
table(eng_contr_PSB=="7") # alta intensidade controle (24/128) = 0.1875
table(eng_ben_PSB=="7") # alta intensidade benef�cio pr�prio partido (54/128) = 0.421875
table(eng_odio_PSB=="7") # alta intensidade �dio (40/128) = 0.3125
table(eng_demint_PSB=="7") # alta intensidade democracia interna (38/128) = 0.296875

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PSB=="7", eng_ben_PSB=="7", alternative = "less", paired = T) #p-value = 7.383e-07

#efeito do engajamento por �dio
t.test(eng_contr_PSB=="7", eng_odio_PSB=="7", alternative = "less", paired = T) #p-value = 0.004458

#efeito do demint
t.test(eng_contr_PSB=="7", eng_demint_PSB=="7", alternative = "less", paired = T) #p-value = 0.00788

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PSB=="7", eng_odio_PSB=="7", alternative = "greater", paired = T) #p-value = 0.00126***


#### PTB ####
dadosPTB <- dados[which(dados$partido=="PTB"),]
eng_contr_PTB <- dadosPTB$hs_controle_cat
eng_contr_PTB <- as.factor(eng_contr_PTB)

eng_ben_PTB <- dadosPTB$hs_benef_cat
eng_ben_PTB <- as.factor(eng_ben_PTB)

eng_odio_PTB <- dadosPTB$hs_odio_cat
eng_odio_PTB <- as.factor(eng_odio_PTB)

eng_demint_PTB <- dadosPTB$hs_demint_cat
eng_demint_PTB <- as.factor(eng_demint_PTB)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PTB) #n=73
table(eng_contr_PTB=="7") # alta intensidade controle (10/73) = 0.1369863
table(eng_ben_PTB=="7") # alta intensidade benef�cio pr�prio partido (25/73) = 0.3424658
table(eng_odio_PTB=="7") # alta intensidade �dio (24/73) = 0.3287671
table(eng_demint_PTB=="7") # alta intensidade democracia interna (20/73) = 0.2739726

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PTB=="7", eng_ben_PTB=="7", alternative = "less", paired = T) #p-value = 0.0006693

#efeito do engajamento por �dio
t.test(eng_contr_PTB=="7", eng_odio_PTB=="7", alternative = "less", paired = T) #p-value = 0.0006622

#efeito do demint
t.test(eng_contr_PTB=="7", eng_demint_PTB=="7", alternative = "less", paired = T) #p-value = 0.01215

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PTB=="7", eng_odio_PTB=="7", alternative = "greater", paired = T) #p-value = 0.3827 (n est sign)


#### PCdoB ####
dadosPCdoB <- dados[which(dados$partido=="PCdoB"),]
eng_contr_PCdoB <- dadosPCdoB$hs_controle_cat
eng_contr_PCdoB <- as.factor(eng_contr_PCdoB)

eng_ben_PCdoB <- dadosPCdoB$hs_benef_cat
eng_ben_PCdoB <- as.factor(eng_ben_PCdoB)

eng_odio_PCdoB <- dadosPCdoB$hs_odio_cat
eng_odio_PCdoB <- as.factor(eng_odio_PCdoB)

eng_demint_PCdoB <- dadosPCdoB$hs_demint_cat
eng_demint_PCdoB <- as.factor(eng_demint_PCdoB)


#descri��o do engajamento de alta intensidade nos cen�rios controle e benef�cio, �dio e demint
length(eng_contr_PCdoB) #n=73
table(eng_contr_PCdoB=="7") # alta intensidade controle (17/73) = 0.2328767
table(eng_ben_PCdoB=="7") # alta intensidade benef�cio pr�prio partido (37/73) = 0.5068493
table(eng_odio_PCdoB=="7") # alta intensidade �dio (28/73) = 0.3835616
table(eng_demint_PCdoB=="7") # alta intensidade democracia interna (23/73) = 0.3150685

#efeito do engajamento em benef�cio do pr�prio partido
t.test(eng_contr_PCdoB=="7", eng_ben_PCdoB=="7", alternative = "less", paired = T) #p-value = 3.006e-06

#efeito do engajamento por �dio
t.test(eng_contr_PCdoB=="7", eng_odio_PCdoB=="7", alternative = "less", paired = T) #p-value = 0.00338

#efeito do demint
t.test(eng_contr_PCdoB=="7", eng_demint_PCdoB=="7", alternative = "less", paired = T) #p-value = 0.05466

#comparando efeito do engajamento em benef�cio do pr�prio partido e o engajamento por �dio 
t.test(eng_ben_PCdoB=="7", eng_odio_PCdoB=="7", alternative = "greater", paired = T) #p-value = 0.005772***



#### Rejei��o ####
rej_algum <- dados$rej_algum
table(rej_algum)
#1195/1917 (0.6233698 disseram n gostar de pelo menos um partido)


rej_PSL <- dados$rej_PSL
table(rej_PSL)
#482/1917 (0.2514345)

rej_PT <- dados$rej_PT
table(rej_PT)
#653/1917 (0.3406364)

rej_MDB <- dados$rej_MDB
table(rej_MDB)
#396/1917 (0.2065728)

rej_DEM <- dados$rej_DEM
table(rej_DEM)
#399/1917 (0.2081377)

rej_PSDB <- dados$rej_PSDB
table(rej_PSDB)
#448/1917 (0.2336985)

rej_PDT <- dados$rej_PDT
table(rej_PDT)
#159/1917 (0.0829421)

rej_PSB <- dados$rej_PSB
table(rej_PSB)
#147/1917 (0.07668232)

rej_PSOL <- dados$rej_PSOL
table(rej_PSOL)
#524/1917 (0.2733438)

rej_PP <- dados$rej_PP
table(rej_PP)
#240/1917 (0.1251956)

rej_PCdoB <- dados$rej_PCdoB
table(rej_PCdoB)
#475/1917 (0.247783)

rej_outros <- dados$rej_outros
table(rej_outros)
#213/1917 (0.1111111)
